Genetic Diversity Analysis of Soybean Collection Using Simple Sequence Repeat Markers.

Soybean [Glycine max (L.) Merr.] is a nutrient-rich crop that offers a sustainable source of dietary protein and edible oil. Determining the level of genetic diversity and relationships between various genetic resources involved in breeding programs is very important in crop improvement strategies. This study evaluated 100 soybean accessions with diverse origins for 10 important agronomic traits, including plant height (PH), an important plant adaptation-related trait impacting yield, in conditions in southeastern Kazakhstan for 2 years. The comparison of different groups of PH (tall, middle, and short) using a t-test suggested that the group of plants with the tallest PH provided a higher yield (p < 0.001) in relatively dry field conditions. The genetic diversity of the accessions was estimated using 25 simple sequence repeat (SSR) markers previously known to be associated with plant height. The results showed a significant variation among different groups of origin for all measured agronomic traits, as well as high genetic diversity, with the PIC (polymorphism information content) varying from 0.140 to 0.732, with an average of 0.524. Nei's diversity index ranged between 0.152 and 0.747, with an average of 0.526. The principal coordinate analysis (PCoA) of the studied soybean collection showed that Kazakhstan accessions were genetically distant from European, East Asian, and North American cultivars. Twelve out of twenty-five SSR markers demonstrated significant associations with ten studied agronomic traits, including PH (p < 0.05). Six SSRs with pleiotropic effects for studied traits were selected, and their haplotypes with phenotypic effects were generated for each soybean accession. The obtained results can be used in soybean improvement programs, including molecular-assisted breeding projects.

Identification of Quantitative Trait Loci Associated with Powdery Mildew Resistance in Spring Barley under Conditions of Southeastern Kazakhstan.

Barley (Hordeum vulgare L.) is one of the most produced cereal crops in the world. It has traditionally been used for the production of animal feed and for malting, as well as for human consumption. However, its production is highly affected by biotic stress factors, particularly the fungal pathogen Blumeria graminis (DC.) f. sp. hordei (Bgh), which causes powdery mildew (PM). In this study, a collection of 406 barley accessions from the USA, Kazakhstan, Europe, and Africa were assessed for resistance to PM over a 3-year period in southeastern Kazakhstan. The collection was grown in the field in 2020, 2021, and 2022 and was genotyped using the 9K SNP Illumina chip. A genome-wide association study (GWAS) was conducted to identify the quantitative trait loci (QTLs) associated with PM resistance. As a result, seven QTLs for PM resistance were detected on chromosomes 4H, 5H, and 7H (FDR p-values < 0.05). Genetic positions of two QTLs were similar to those of PM resistance QTLs previously reported in the scientific literature, suggesting that the five remaining QTLs are novel putative genetic factors for the studied trait. Haplotype analysis for seven QTLs revealed three haplotypes which were associated with total PM resistance and one haplotype associated with the high PM severity in the barley collection. Identified QTLs and haplotypes associated with the PM resistance of barley may be used for further analysis, trait pyramiding, and marker-assisted selection.

Genome-Wide Association Study of Leaf Rust and Stem Rust Seedling and Adult Resistances in Tetraploid Wheat Accessions Harvested in Kazakhstan.

Leaf rust (LR) and stem rust (SR) are diseases increasingly impacting wheat production worldwide. Fungal pathogens producing rust diseases in wheat may cause yield losses of up to 50−60%. One of the most effective methods for preventing such losses is the development of resistant cultivars with high yield potential. This goal can be achieved through complex breeding studies, including the identification of key genetic factors controlling rust disease resistance. The objective of this study was to identify sources of tetraploid wheat resistance to LR and SR races, both at the seedling growth stage in the greenhouse and at the adult plant stage in field experiments, under the conditions of the North Kazakhstan region. A panel consisting of 193 tetraploid wheat accessions was used in a genome-wide association study (GWAS) for the identification of quantitative trait loci (QTLs) associated with LR and SR resistance, using 16,425 polymorphic single-nucleotide polymorphism (SNP) markers in the seedling and adult stages of plant development. The investigated panel consisted of seven tetraploid subspecies (Triticum turgidum ssp. durum, ssp. turanicum, ssp. turgidum, ssp. polonicum, ssp. carthlicum, ssp. dicoccum, and ssp. dicoccoides). The GWAS, based on the phenotypic evaluation of the tetraploid collection's reaction to the two rust species at the seedling (in the greenhouse) and adult (in the field) stages, revealed 38 QTLs (p < 0.001), comprising 17 for LR resistance and 21 for SR resistance. Ten QTLs were associated with the reaction to LR at the seedling stage, while six QTLs were at the adult plant stage and one QTL was at both the seedling and adult stages. Eleven QTLs were associated with SR response at the seedling stage, while nine QTLs were at the adult plant stage and one QTL was at both the seedling and adult stages. A comparison of these results with previous LR and SR studies indicated that 11 of the 38 QTLs are presumably novel loci. The QTLs identified in this work can potentially be used for marker-assisted selection of tetraploid and hexaploid wheat for the breeding of new LR- and SR-resistant cultivars.

Population Structure and Genetic Diversity of Two-Rowed Barley Accessions from Kazakhstan Based on SNP Genotyping Data.

The genetic relationship and population structure of two-rowed barley accessions from Kazakhstan were assessed using single-nucleotide polymorphism (SNP) markers. Two different approaches were employed in the analysis: (1) the accessions from Kazakhstan were compared with barley samples from six different regions around the world using 1955 polymorphic SNPs, and (2) 94 accessions collected from six breeding programs from Kazakhstan were studied using 5636 polymorphic SNPs using a 9K Illumina Infinium assay. In the first approach, the neighbor-joining tree showed that the majority of the accessions from Kazakhstan were grouped in a separate subcluster with a common ancestral node; there was a sister subcluster that comprised mainly barley samples that originated in Europe. The Pearson's correlation analysis suggested that Kazakh accessions were genetically close to samples from Africa and Europe. In the second approach, the application of the STRUCTURE package using 5636 polymorphic SNPs suggested that Kazakh barley samples consisted of five subclusters in three major clusters. The principal coordinate analysis plot showed that, among six breeding origins in Kazakhstan, the Krasnovodopad (KV) and Karaganda (KA) samples were the most distant groups. The assessment of the pedigrees in the KV and KA samples showed that the hybridization schemes in these breeding stations heavily used accessions from Ethiopia and Ukraine, respectively. The comparative analysis of the KV and KA samples allowed us to identify 214 SNPs with opposite allele frequencies that were tightly linked to 60 genes/gene blocks associated with plant adaptation traits, such as the heading date and plant height. The identified SNP markers can be efficiently used in studies of barley adaptation and deployed in breeding projects to develop new competitive cultivars.

Identification of Quantitative Trait Loci for Leaf Rust and Stem Rust Seedling Resistance in Bread Wheat Using a Genome-Wide Association Study.

In recent years, leaf rust (LR) and stem rust (SR) have become a serious threat to bread wheat production in Kazakhstan. Most local cultivars are susceptible to these rusts, which has affected their yield and quality. The development of new cultivars with high productivity and LR and SR disease resistance, including using marker-assisted selection, is becoming an important priority in local breeding projects. Therefore, the search for key genetic factors controlling resistance in all plant stages, including the seedling stage, is of great significance. In this work, we applied a genome-wide association study (GWAS) approach using 212 local bread wheat accessions that were phenotyped for resistance to specific races of Puccinia triticina Eriks. (Pt) and Puccinia graminis f. sp. tritici (Pgt) at the seedling stages. The collection was genotyped using a 20 K Illumina iSelect SNP assay, and 11,150 polymorphic SNP markers were selected for the association mapping. Using a mixed linear model, we identified 11 quantitative trait loci (QTLs) for five out of six specific races of Pt and Pgt. The comparison of the results from this GWAS with those from previously published work showed that nine out of eleven QTLs for LR and SR resistance had been previously reported in a GWAS study at the adult plant stages of wheat growth. Therefore, it was assumed that these nine common identified QTLs were effective for all-stage resistance to LR and SR, and the two other QTLs appear to be novel QTLs. In addition, five out of these nine QTLs that had been identified earlier were found to be associated with yield components, suggesting that they may directly influence the field performance of bread wheat. The identified QTLs, including novel QTLs found in this study, may play an essential role in the breeding process for improving wheat resistance to LR and SR.

Genome-wide association mapping for resistance to leaf, stem, and yellow rusts of common wheat under field conditions of South Kazakhstan.

Common or bread wheat (Triticum aestivum L.) is the most important cereal crop in the world, including Kazakhstan, where it is a major agricultural commodity. Fungal pathogens producing leaf, stem, and yellow (stripe) rusts of wheat may cause yield losses of up to 50-60%. One of the most effective methods for preventing these losses is to develop resistant cultivars with high yield potential. This goal can be achieved using complex breeding studies, including the identification of key genetic factors controlling rust disease resistance. In this study, a panel consisting of 215 common wheat cultivars and breeding lines from Kazakhstan, Russia, Europe, USA, Canada, Mexico, and Australia, with a wide range of resistance to leaf rust (LR), stem rust (SR), and yellow rust (YR) diseases, was analyzed under field conditions in Southern Kazakhstan. The collection was genotyped using the 20K Illumina iSelect DNA array, where 11,510 informative single-nucleotide polymorphism markers were selected for further genome-wide association study (GWAS). Evaluation of the phenotypic diversity over 2 years showed a mostly mixed reaction to LR, mixed reaction/moderate susceptibility to SR, and moderate resistance to YR among wheat accessions from Kazakhstan. GWAS revealed 45 marker-trait associations (MTAs), including 23 for LR, 14 for SR, and eight for YR resistances. Three MTAs for LR resistance and one for SR resistance appeared to be novel. The MTAs identified in this work can be used for marker-assisted selection of common wheat in Kazakhstan in breeding new cultivars resistant to LR, SR, and YR diseases. These findings can be helpful for pyramiding genes with favorable alleles in promising cultivars and lines.

Association mapping for agronomic traits in six-rowed spring barley from the USA harvested in Kazakhstan.

In barley, six-rowed barley is advantageous over two-rowed barley for feed due to the larger number of seeds per spike and the higher seed protein content. The growth of six-rowed barley is potentially important for breeding in agriculturally oriented countries, such as Kazakhstan. Nevertheless, until recently, very little attention was given to six-rowed barley in breeding projects in Kazakhstan, one of the largest countries in the world. In this study, phenotyping and single nucleotide polymorphism (SNP) genotyping data were generated from 275 accessions originating from six different breeding organizations in the USA as well as 9 accessions from Kazakhstan in field trials at six breeding institutions. The USA six-rowed barley was tested in comparison to local accessions over three years (2009-2011) based on analyses of key agronomic traits. It was determined that the average yield in the USA accessions in comparison to local lines showed heavier yield in all six tested sites. Principal Coordinate Analysis based on 1618 polymorphic SNP markers separated Kazakh lines from six USA barley origin groups based on PC1 (77.9%), and Montana lines from the remaining five USA groups based on PC2 (15.1%). A genome-wide association study based on eighteen field trials allowed the identification of 47 stable marker-trait associations (MTA) for ten agronomic traits, including key yield related characters such as yield per square meter, thousand grain weight, number of kernels per spike, and productive tillers. The comparison of chromosomal positions of identified MTA with positions of known genes and quantitative trait loci suggests that 25 out of those 47 MTAs are presumably novel. The analysis of 42 SNPs associated with 47 MTAs in the Ensemble genome annotation system (http://ensemblgenomes.org) suggested that 40 SNPs were in genic positions of the genome, as their sequences successfully aligned with corresponding Gen ID.

Marker-trait associations in two-rowed spring barley accessions from Kazakhstan and the USA.

In this study, phenotyping and single nucleotide polymorphism (SNP) genotyping data of 272 accessions of two-rowed spring barley from the USA along with 94 accessions from Kazakhstan were assessed in field trials at six breeding organizations in Kazakhstan to evaluate the performance of the USA samples over three years (2009-2011). The average grain yield over the six locations was not significantly higher in Kazakh accessions in comparison to the USA samples. Twenty four samples from Montana, Washington, the USDA station in Aberdeen Idaho, and the Anheuser-Busch breeding programs showed heavier average yield than the local standard cultivar "Ubagan". Principal Coordinate analysis based on two sets of SNP data suggested that Kazakh accessions were closest to the USA accessions among eight groups of samples from different parts of the World, and within five US barley origin groups the samples from Montana and Washington perfectly matched six groups of Kazakh breeding origins. A genome-wide association study (GWAS) using data from eighteen field trials allowed the identification of ninety one marker-trait associations (MTA) in two or more environments for nine traits, including key characters such as heading time (HT), number of kernels per spike (NKS), and thousand grain weight (TGW). Our GWAS allowed the identification of eight MTA for HT and NKS, and sixteen MTA for TGW, when those MTA were linked to mapped SNPs. Based on comparisons of chromosomal positions of MTA identified in this study, and positions of known genes and quantitative trait loci for HT, NKS and TGW, it was suggested that MTA for HT on chromosome 2H (at 158.2 cM, 11_21414), MTA for NKS on 5H (at 118.6 cM, 11_20298), and two MTA for TGW on chromosome 4H (at 94.7 cM, 12_30718, and at 129.3 cM, 11_20013) were potentially new associations in barley. GWAS suggested that six MTA for HT, including two on chromosome 1H, two on chromosome 3H, and one each on chromosomes 4H and 6H, had useful pleiotropic effects for improving barley spike traits.

Morphological description and DNA barcoding study of sand rice (Agriophyllum squarrosum, Chenopodiaceae) collected in Kazakhstan.

BACKGROUND: Sand rice (Agriophyllum squarrosum (L.) Moq.) is an annual shrub-like plant adapted to the mobile sand dunes in desert and semi-desert regions of Asia. It has a balanced nutrient composition with relatively high concentration of lipids and proteins, which results in its nutrition being similar to legumes. Sand rice's proteins contain the full range of essential amino acids. However, calories content is more similar to wheat. These features together with desert stress resistance make sand rice a potential food crop resilient to ongoing climate change. It is also an important fodder crop (on young stages of growth) for cattle in arid regions of Kazakhstan. In our work, sand rice samples were collected from two distant regions of Kazakhstan as a part of the nation-wide project to determine genetic variation of the native flora. RESULTS: Samples were collected in western and southeastern parts of Kazakhstan separated by distances of up to 1300 km. Sequences of the nuclear ribosomal DNA ITS1-5.8S-ITS2 region and the chloroplast matK gene confirmed the identity of species defined by morphological traits. Comparison with GenBank sequences revealed polymorphic sequence positions among Kazakh populations and GenBank references, and suggested a distinction among local populations of sand rice. The phylogenetic analysis of nucleotide sequences showed a clear partition of A. squarrosum (L.) Moq. from Agriophyllum minus Fisch. & C.A. Mey, which grows in the same sand dunes environment. CONCLUSIONS: DNA barcoding analyses of ITS and matK sequences showed a segregation of A. squarrosum from A. minus into separate clades in Maximum-Likelhood dendrograms. ITS analysis can be successfully used to characterize A. squarrosum populations growing quite distant from each other. The data obtained in this work provide the basis for further investigations on A. squarrosum population structure and may play a role in the screening of sand rice plants growing in desert and semi-desert environments of Central Asia and China.

Taxonomic assessment of Allium species from Kazakhstan based on ITS and matK markers.

BACKGROUND: As part of nation-wide project to infer the genetic variation of the native flora in Kazakhstan, a study was attempted to assess phylogenetic relationships of endemic and rare Allium species. In total, 20 Allium species were collected in field trips in five different regions of Kazakhstan during 2015-2016. Most species (9) were collected in the southern part of the country along of Karatau mountains, followed by Altai mountains (5) in eastern Kazakhstan. The ITS and matK DNA regions were applied in order to assess the taxonomic relationships among species. The major goal of the study was to assess the taxonomic position of five endemic and rare species from Allium subgenus Reticulatobulbosa collected in Karatau mountains of Southern Kazakhstan. RESULTS: The 20 collected Allium species were assessed using morphological traits and a DNA barcoding approach. The morphological analyses of four different species in subgenus Reticulatobulbosa inferred similarities of A. inconspicuum and A. barszchewskii (both from section Companulata) that were separated from A. oreoscordum and A. oreoprasoides (section Nigrimontana) by several traits, including form of bulbs and leaves, presence of bracts, shape of perianth lobes and style. The Neighbor-Joining method was applied to generate ITS and matK phylogenetic trees for two groups of populations: 1) 20 Allium species collected within the project, and 2) 50 Allium worldwide species. CONCLUSIONS: The analyses of nucleotide sequences of ITS and matK robustly confirmed the monophyletic origin of the Allium species. The variability in 20 local Allium species in ITS was 6.6 higher than in matK, therefore the topology of the ITS tree was better resolved. The taxonomy of Allium species largely coincided with a recent classification of this genus. Analyses of both ITS and matK suggest that A. oreoscordum is genetically close to A. oreoprasoides in section Nigrimontana of subgenus Reticulatobulbosa. This result was also confirmed using morphological description of individual plants of four species in subgenus Reticulatobulbosa. The study is another contribution to taxonomy clarification in Allium.